Exact vortex nucleation and cooperative vortex tunneling in dilute Bose-Einstein 

Condensates 
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With the imminent advent of mesoscopic rotating BECs in the lowest Landau level (LLL) regime, 
we explore LLL vortex nucleation. An exact many-body analysis is presented in a weakly elliptical 
trap for up to 400 particles. Striking non-mean field features are exposed at filling factors ^ 1. 
E. g. near the critical rotation frequency pairs of energy levels approach each other with exponential 
accuracy. A physical interpretation is provided by requantising a mean field (MF) theory, where 1/N 
plays the role of Planck's constant, revealing two vortices cooperatively tunneling between classically 
degenerate energy minima. The tunnel splitting variation is described in terms of frequency, particle 
number and ellipticity. 

PACS numbers: 03.75.Hh, 03.75.Lm 
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The physics of vortices in slowly rotating degenerate 
gases [l| has reached the level of maturity where it is 
now used as a tool to study other phenomena, such as 
polarised fermi gases 0- However achieving rapid rota- 
tion - to explore thoroughly the MF quantum Hall (QH) 
regime[l, U 0, 0] in the lowest Landau level (LLL) 
and to reach corr elated QH states 0, i, % - remains a 
challenge. 

A promising approach to accessing the QH regime is 
to have very dilute BECs, perhaps constructed by slicing 
up a condensate with an optical lattice[l^. In this Letter 
we show that even well away from the correlated regime 
there are pronounced quantum effects which become in- 
creasingly striking as the particle number decreases. We 
will show that the exact many-body ground states may 
be interpreted as exhibiting vortex tunneling leading to 
superpositions of mean-field states with vortices residing 
at different locations. This mesoscopic limit is consistent 
with the thrust of experimental effort in the near future. 
(In terms oi v = N/N^, where N is the number of parti- 
cles, and the number of vortices, = 1/2 corresponds 
to the Laughlin state, and we will study 10 < < 400.) 



Jf, for a cold gas of N particles residing in a plane: 



Vortex nucleation [111 has been studied in the Thomas- 
Fermi regime, both experimentally [l^ 13 . 14 1 and 
theoretically [Ta. The conclusion is that under adi- 
abatic ramping of the rotation frequency [l^, the 
process is determined by an hydrodynamic instability. 
Under those conditions, the thermodynamic instability 
to vortex entry is apparently unobservable, occurring at 
lower rotation frequencies. 

It is known [TtI that in a BEC in the LLL in an ax- 
isymmetric(AS) trap that there is a first-order thermo- 
dynamic instability to vortex entry (with no hydrody- 
namic instability needed). In this Letter, we will show 
that the situation is very different in a non-AS trap. The 
equilibrium of vortices in a non-AS trap has already been 
analysed at a MF level[liJi3,[2i within the LLl[2i| and 
at the Bogohubov levelj22|). 

Our starting point is the standard model Hamiltonian, 
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Units of length, aj_, and energy, hu!j_, arc those provided 
by the harmonic trap; angular momenta, are scaled 
by H. There are two remaining dimensionless parameters. 
Firstly, fi, is the angular velocity of the rotating frame 
divided by the natural frequency of the harmonic trap. 
Secondly the couphng constant, 77 = ATra/a±. We assume 
the particles reside in the LLL of this AS trap 01 ■ The 
LLL single-particle basis utilises a complex description 
of the particle positions in the plane, z = x + iy, and 
is defined by the set {(z'^/a/tttoI) e"'^' /^}, where m ~ 

o,i,---. 

In the AS case the many-body eigenfunctions, 
'>pL{{zk}), (labeled by the total angular momentum, L) 
are known analytically for i < 0, [13, HI] ■ They are 



1pL{{zk})0C {Z^^-Zc){z^^ 
l<ii---iL<N 



Zc^ ' ' ' {Zij^ Zc) 



where Zc = ^/N'^^^^Zj is the centre of mass co- 
ordinate. The energies, E^, are also known: for < L < 
N{L 7^ 1),(L = 1 is special as it corresponds to centre of 
mass motion and the energy is = 1 — fi). 



El = N + L + \7^N{N -\L-1)-VlL. 
= N+ \riN{N - 1) + (r2(°) - VL)L. 



(1) 
(2) 



At the critical frequency, fic"^ = 1—^, all the eigenstates 
for 2 < i < A^, and L = 0), are degenerate. The entry is 
abrupt, but degenerate. 

In this Letter we determine how that degeneracy is 
lifted in a weakly elliptical trap with striking physical 
consequences. The elliptical perturbation has the form: 
\e'K2 = \ f-Y^n=i i^n + ^n^)- Wc wiU assume that e<.rj, 
to allow a description in terms of the AS states. 
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FIG. 1: The low-lying energy gaps, AE„, measured from the 
ground state for A*' = 400 when the trap is rotated in the 
vortex formation region near f2 ~ 0. Inset; detail of the low- 
lying energy gaps near Qc, when the degeneracies lift. 

To expose clearly the behaviour in the vicinity of vor- 
tex nucleation, it is convenient to use equation ^ . Then, 
changing the zero of energy to absorb the term indepen- 
dent of L, we find the complete rescaled Hamiltonian 
(choosing e > as the unit of energy) 

where Cl = {fl — fli^^)/e. This rescaling stretches out 
the nucleation region where our approximation is valid 
(e <§; rj) which is convenient numerically. When the en- 
ergy levels are calculated by exact diagonalisation (ED) 
for TV = 400, we find the results portrayed in FiglT] (the 
energies are measure relative to the ground-state energy 
and we scale rj = r]o/N to provide a sensible thermody- 
namic limit). We restrict ourselves to even N, as for odd 
N there is a trivial residual first order transition (due to 
the necessary change of parity of the ground state) whose 
magnitude diminishes as ^ oo. 

To understand the significance of the results, it is 
useful to consider what the equivalent diagram would 
have shown in the AS case. There, from equation ([1]), 
we should measure the energy from the ground state 
Eg ~ El=o for Q < Oc, and from Eg = El=n for 
r2 > r^c, i.e. the appropriate ground states (within the 
basis set in the latter case at higher J7). Then we find: 

^ ^ { {Q - n^){N - L) : fl~flc>0 

This would also give a "V" -shape, however there are sub- 
stantial differences. Firstly, in the AS case, the apex of 
the "V" is at r2 — J7c = 0, whereas the elliptical case is 
displaced by roughly e/4. Secondly, in the AS case the 
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FIG. 2; The tunneling coefficient a, fitted to iV' e""^ from 
160-450 particles. On the left hand side of the peak (3 ~ 0.53 
but the right hand side does not scale simply. 

gradients increase linearly with level index. However in 
FiglT] the energy levels form doublets (which are of oppo- 
site parity) as they approach the critical point, and the 
doublets themselves do not have simply related gradients 
from doublet to doublet. Of course, the expanded scale 
of the figure (due to being scaled by e) emphasises this 
and at a sufficient distance from the critical point, the 
AS gradients must be obtained. 

Focussing on the minimum (as a function of CI) gap, 
we have performed a finite size scaling fit of the gap be- 
tween the ground state and first excited state, in terms 
of the dependence on the number of particles, N, to 
AE = N^e~°'^ . This form is consistent with a tunneling 
process (which the doublet structure suggests) occurring 
prior to the vortex being in the centre of the BEG. Data 
for Ai? > 10"^'' has been kept for a range of N and a fit 
found for a, as seen in FigO 

The most significant feature, is a peak in Q;(f2), whose 
height appears constant with increasing TV. Having found 
the behaviour of a, we find that the choice of the ex- 
ponent (3 K, 0.53 collapses the data for all N onto one 
curve. On the right hand side, there appears to be no 
single choice for /3, to collapse the data onto one curve. 

The above analysis suggests, via the exponential de- 
pendence on TV and the pairs of opposite parity eigen- 
states, that tunneling is involved in the nature of the 
states near ilc- By analogy with a particle in a double 
well, the states linked by tunneling can be exposed by 
examining the states V'i = l/'\/2(V'o ± "01); which would 
correspond to states, in the double well case, where the 
particle is localised in one or the other well. The single 
particle densities, p±, associated with these states, reveal 
a depression of the density on one side (which is deter- 
mined by which of the states is being considered) of the 
semi-major axis. This depression in density moves to- 
wards the centre of the trap as O ^ f2c • Further analysis 
of Fig.[l]indicates that gaps between higher pairs of oppo- 
site parity levels also vanish asymptotically as TV ^ oo. 

To establish a physical interpretation of these results, 
we extend the variational LLL MF Lagrangian[2^. to the 
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elliptical case. We work with two vortices, as this is the 
simplest system with enough parameters to encapsulate 
the qualitative features of the exact solution. Following 
2J| we use hydrodynamic variables, {pm, 4'm}, such that 
the variational state ipMiz,z*) = X]rn=i ^™^™®~'^ '^^ 
and ttm = \/ PTn/(7rTO!)e"^™. The iwo- vortex Lagrangian 
corresponds to maximum angular momentum M ~ 2. In 
units of e, the two vortex Lagrangian becomes (in the 
first line in unsealed form): 



= hN[p^cb_+p+(j)+ 

21 
e 



r(!p 



r - O ) (1 - p_) + cr(/o+, p_) cos(2</._) 

1 
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2(l-2p+)a(p+,p_)cos0+). (3) 



where we have: transformed to variables p+ = ^(Po+P2), 
0+ = (00 + 02), p- = (po - P2) and 0_ = 1(00 - 02); 
defined ct(/9+,P-) = (l/(2^/2))(4p2_ - p2^)i/2 and T = 
77o/4e7r^-\/27r; and used normalisation X]i=o Pi = 1, to 
eliminate pi. (In the vicinity of vortex nucleation pi ~ 
1 and po and p2 are correspondingly small.) Wc have 
picked the arbitrary phase such that 0i = 0. 

At frequencies, 17 < f^c, consideration of the La- 
grangian indicates two degenerate energy minima. By in- 
spection for the phase variables, 0±, we see that 0+ = tt 
and 0_ = ±7r/2. Minimising the Hamiltonian numeri- 
cally with respect to the p± gives the equilibrium values 
p± for a given set of interaction and rotation parame- 
ters. In the original variables, the minima have associ- 
ated phase variables 0o = 0, 02 = tt, and 0o = tt, 02 = 0. 
This implies that both minima correspond to two vortices 
on the X-axis (one cither side of the origin, not necessar- 
ily at equal distance) and are mirror refiections of each 
other. Beyond MF one might expect that the vortices tun- 
nel between these configurations if they are not identical. 

Noting that the potential involving cos0_ is the only 
term not multiplied by the large parameter, F, we Taylor 
expand all the other terms to harmonic level around the 
equilibria. We also evaluate cr, at the equilibrium values 
of p± leading to a potential ~ ea{p+, p^ ) cos cj) 

We follow the standard procedure e.g. Ref. |26|, to 
examine quantum effects beyond MF theory, and re- 
quantise the Lagrangian, which now only involves the 
— variables, and start from Eqn. ([3]). The classical con- 
jugate variables to 0± are p± = Nhp±. Thus the com- 
mutator for canonically conjugate variables: 



i>^,Nhp^] 



ih 



=> [4>-,p-\ 



Hence the usual quantisation procedure implies: /5_ 
—\N~^d/d(j)- This is clearly unaltered under the dis- 
placement of p- = p- + Sp- . We note that the quantum 
effects will vanish in the limit N 00, i.e. the thermo- 
dynamic limit, consistent with the exact results described 
above. 



Making the operator replacements, leads to 
^ '^'--N^£~aip+,p-)cosi2(b-)) 



2M* 



1'(0-) = 

(4) 

where we have defined: 1/M* = d^T/ddp^\sp_=o and 
£ ~ [E ~ Vb)/iV, with Vb being the zcroth order terms 
from the T expansion. This is a Mathieu equation. 



[dxx + b- 2gcos(2a:;)] y{x) = 



(5) 



and inspection shows that h = 2N'^M*£ > and q = 
N'^M*a{p+, P-) > 0. In the limiting case of g = 0, the 
Mathieu functions are simply cos(a/60_) and sm{Vb(j)-) , 
so the AS result is recovered. 

From Eqn. (20.2.31) [111, the level splitting (between 
the first excited and ground states) of Eqn. ((51) is 
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implying 



iV4e' 



where the tunnelling coefficient 



that A£' 

a2v = A^M*a{p+,p-). 

Finally reinstating the unsealed 77 we can estimate the 
tunnel splitting using say = 10, e = 10^'^ and ^ ~ 1 
(i.e. a chip trap). We then find E ~ O.S^kjQ, indicating 
the splitting should be observable under those assump- 
tions. 

These results reproduce the exponential suppression 
of the tunneling splitting of the eigenvalues found in the 
ED, indicating we have correctly identified the tunnel- 
ing entities. Qualitatively the behavior of the tunneling 
coefficient, a, with rotation frequency (shown in Fig. ^ 
follows that of the exact result, having a maximum in 
Q!(0). We have scaled the MF frequencies, O, by 1/4 
to compare the peaks in a. This reflects ~ e 

as against ^7°^^^=* = g/4. In addition the exponent of 
the pre-exponential factor, A^^, has /3™'=^"f = 1/2 (scaled 
variables) almost equal to the ED result, /3°'^'"^* ~ 0.53, 
for Vl to the left of the peak. We found that adding ad- 
ditional (up to eight) vortices improves the agreement 
between the approximate calculation and the exact one, 
at the cost of diminishing the clarity of interpretation 
(the motivation for the requantised calculation). Finite 
temperature will blur the peak, but leave it visible under 
typical experimental conditions. 

Working with two vortices the tunneling path may be 
interpreted in terms of vortex positions. They are: 




i± 1 



2\/2pop2 
Pi 



l/2\ 



Varying 02 along [0,7r] results in both ^± having semi- 
circular trajectories, as shown in Fig. ID Note, that there 
is always a second anticlockwise trajectory correspond- 
ing to 02 from [tt, 27r], and this will result in the major 
semi-circle being above axis. One might have expected 
that the tunneling matrix clement for a vortex would 
increase as the initial and final positions of the vortex 
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FIG. 3: The tunneling parameter from ED (line) and for 
two-vortex MF theory with F = 50, 100, 200 (triangle,star, 
square). Our approximations become more valid for increas- 
ing F. The MF frequency has been scaled to enable qualitative 
comparison of the peaks. 
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FIG. 4: Semicircular trajectories of the vortex positions, en- 
suring the separation between vortices is maintained. 



approached each other (i.e. the inner circle in Fig. 4). 
The more surprising aspect is that the vortex at large 
distances from the centre of the trap does not substan- 
tially decrease the matrix element. This is due to the 
vortex states (labeled by position) becoming increasingly 
non-orthogonal as the vortices move away from the cen- 
tre of the trap. The non-orthogonality arises because the 
Gaussian weight strongly sup presses the region where the 
states are most different [27[ . Hence tunneling is least ef- 
fective when both vortices are at the same, intermediate, 
distance from the trap centre. Exploration of experimen- 
tal observation schemes is underway, including oscillation 
of the trap along the semi-minor axis and noise spectra 
which will reveal the density-density correlation function. 

In conclusion, we have shown, by exact solution, that 
mesoscopic rotating BECs will show pronounced devia- 
tions from MF theory at low angular velocities. This may 



be interpreted, using a requantised MF theory, in terms 
of vortex tunneling in the nuclcation process. 
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